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ABSTRACT 

We present an analytical model of the relation between the surface density of gas and star formation 
rate in galaxies and clouds, as a function of the presence of supersonic turbulence and the associated 
structure of the interstellar medium. The model predicts a power-law relation of index 3/2, flattened 
under the effects of stellar feedback at high densities or in very turbulent media, and a break at 
low surface densities when ISM turbulence becomes too weak to induce strong compression. This 
model explains the diversity of star formation laws and thresholds observed in nearby spirals and 
their resolved regions, the Small Magellanic Cloud, high-redshift disks and starbursting mergers, as 
well as Galactic molecular clouds. While other models have proposed interstellar dust content and 
molecule formation to be key ingredients to the observed variations of the star formation efficiency, 
we demonstrate instead that these variations can be explained by interstellar medium turbulence and 
structure in various types of galaxies. 

Subject headings: stars: formation — galaxies: ISM — hydrodynamics — methods: analytical 



1. INTRODUCTION 

Galactic-scale star formation (SF) laws are not univer- 
sal. Galaxy mergers, when they experience a starburst 
phase, convert their gas into stars within a d epletion time 
up to t en times shorter than spiral galaxies dDaddi et al.l 
I2010ri iGenzel et aLllMot ISaintonge et al.ll2012D . Oppo- 
sitel y, SF in dwarf galaxies is less efficient than in spi- 
rals (jLerov et al.l 120081: iBolatto et al.l 1201 On smaller 
scales, molecular clouds may follow similar scaling laws 
but w ith shorter conversion timescales (e.g. lLada et "all 

To explain these environment-dependent differences, 
some models have emphasized the role of interstellar 
medium (ISM) chemistry (metallicity and dust content) 
in triggering gas cooling and molecule formation. The 
transition from long to short gas depletion times may cor- 
respond to the critical column-density nee ded to shield 
gas from the ambient ultraviolet radiati on (lSchavell2004D 
and/ or to efficiently convert Hi into H2 (jKrumholz et al.l 
2009), which strongly depends on the metallicity. This 
is supported by recent observations of the Small Mag- 
ellanic Cloud (SMC), of lower metallicity than spirals, 
where the star formation rate (SFR) surface density 
hardly reaches the standar d regime even at r elatively 
high gas surface densities (Bo latto et al.l 120111 ). How- 
ever, such descriptions do not explain the very long gas 
depletio n timescales in metal- rich and molecule-rich el- 
lipticals (|Saintonge et al.l2012f) . Furthermore, significant 
SF activity can be found in regions wh ere the molec- 
ular fraction is low (jBoissier et al.l [20081. More gener- 
ally, variations in the SF scaling-laws about as strong in 
H 2 as in total gas have been detected, suggesting that 
molecule formation alo ne does not drive these variations 
(jSaintonge et al.ll2012l ). 

Nevertheless, other factors may explain the observed 
scaling laws. The gas surface density observed in an en- 
tire galaxy or a relatively large region is a global property 
that results from a very heterogeneous distribution of the 
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volume density on small scales. The latter varies broadly 
from quasi-empty holes to dense clouds and cores, un- 
der the effect of the ISM turbulence , which is supersonic 
for a la rge fraction of the mass (e.g. lAudit fc Hennebelld 
I2010D . lElmegreenl (|2002D has shown that th e turbulent 
struct ure of the ISM can naturally explain the lKennicuttl 
( 1998) relation for spiral galaxies. 

In this paper, we develop an analytic model describing 
the variations and thresholds in SF laws depending on 
the supersonic nature of the turbulence and the resulting 
structure of the ISM in several types of galaxies. 

2. THEORY 

2.1. Analytical formalism 

We consider a region of surface S and thickness h, 
whose total mass M is distributed according to a mass- 
weighted probability density function (PDF) /. The vol- 
ume of the region can be written as the sum of the vol- 
umes occupied by the gas at all possible densities p, i.e. 



hS = 



Mf(x) 



dx = — 
P 



dx, 



(1) 



10 P P JO 

where x = p/p is a normalization of the local gas volume 
density p (i.e. for scales < 1 pc), related to gas surface 
density of the region £ = M/S through 



dx. 



The SFR surface density reads 

1 f°° Mf(x) 
S 



J SFR 



PSFR dx, 



(2) 



(3) 



/o P 

where psfr is the local SFR volume density. Therefore, 
the general expression of the SFR surface density ifQ 



^SFR 



Jq 00 f( x ) x 1 Psfr da; 



(4) 



Jo ]_ dx 

1 S ee also Padoan & Nordlund 2011 and Hcmicbcllc fc Chabrierl 
I201U for comparable approaches. 
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2.1.1. Log-normal PDF 

ISM mass is supersonically tur- 

e.g - 

Vazquez-Semadenil 119941 : iNordlund fe Padoanl 119991: 
Wada fe Normanll2001ft : 



Most of the 
bulent, which generate s a log-normal PDF 



fa(x) = 



1 



",\/2<J 2 TT 



exp 



In (s) - %- 
2^2 



(5) 



whose dimensionless width a is related to the Mach num- 
ber M. (and thus to the gas velocity dispersion) through 
a 1 rj In (1 + 3M 2 /4). This was primarily found for 
isothermal gas, but models of non-isothermal ISM in- 
cluding gravity and stellar feedback have shown that de- 
viations from the log-normal functio nal form are negligi- 
ble, at least in isolated disk galaxies (Tas ker fe Tanll2009l : 
iBournaud et~al1l2011bO . 

2.1.2. Dynamical star formation 

At small-scale, cold enough gas becomes supersonically 
turbulent and hosts the shocks which trigger the process 
of star formation. This typically happens above a cer- 
tain volume density po: m nearby spirals for instance, 
the velocity dispersion of 6 — 10 km s _1 requires cooling 
below ~ 10 4 K for ISM turbulence to become supersonic. 
This only becomes possible above po « 10 cm -3 , at so- 
lar metallicty, accordin g to calculations at galactic scale 
(jBournaud et al.H2010D , and detailed ISM models includ- 
ing developed turbulence and radiative transfer which 
show M (p/W cm" 3 ) 1/2 (lAudit fe Hennebelld I20T0L 
Fig. 4 and 9). Above this density threshold, a fraction of 
the gas becomes gravitationally unstable and collapses, 
converting a constant fractional e = 0.01 of its mass into 
stars per free-fall time tg, as supported by observations 
(|Elmegreenll2002D : 





/'Sl-'R { P_ 

C tg 



2.1.3. 



if p < Po 



32G 
"3^ 



P 



• 3/2 else 



(6) 



Regulation by stellar feedback 

In dense cores, stellar feedback limits the conversion of 
the gas into stars, by heating, ionizing and even eject- 
ing the gaseous left-overs. This gas becomes available 
again for SF after a time t s . No more than e s = 30% of 
the gas mass can be consumed for SF per time -scale t s 
(|Bontemps et al.|[l99l [Matzner fe McKeel |2000T) . This 
translates into a saturation in the local SF law: 







/'sm; \ min ( e ,/^ p 3/2 e P 

J7T t s 



if P < Po 
else ' 



The (re) formation of star forming cloud s is generally trig- 
gered by galactic- scale processes (e.g. iDobbs fe Pringld 



20091 for spirals iTevssier et all |2010| for mergers and 
Bournaud et al.l 120071 at high redshift). therefore t s 



Observations suggest e as 0.01, inde pendently of the local 
density and scale (Krumholz & Tan 2007). Wc adopt this value 
throughout the paper. 



100 Myr over a broad range of galaxy masses. For in- 
stance, this represents the interval between the compres- 
sion by two spiral arms, or for giant clumps to collapse 
in high-redshift disks. In nuclear starbursts however, the 
dynamical timescale is much shorter (~ 10 Myr) but 
strong stellar feedback from OB-type stars takes over 
and also regulates SF , by limiting the co nversion of gas 
into stars to e s sa 0.3 ((Murray et al.ll2010l ) over the dura- 
tion of the starburst ev ent, i.e. t s = 100 Myr here again 
(|Di Matteo et al.ll2008f ). 

2.2. Results 

Without feedback (equation [6|) and for a log-normal 
PDF, we obtair0 



J SFR 




S 3/2 erfc 



In 



( Po h \ 

v s ; 



(8) 

As illustrated in Fig. [TJabc, this corresponds to a 3/2- 
index power law which, for any non-zero threshold pq, 
falls off at low surface densities. Hence, the combination 
of standard functional forms for the density PDF of the 
ISM and the local SFR naturally results in a threshold 
(or "break") in Ssfr, the shape and position of which 
depend on A4, h and pq. 

As shown in Fig.[T]d, the regulation by feedback (equa- 
tion [7]) induces a shallower power law (of index unity) at 
high surface densities: 



Ssfr — °a,o S 3//2 + c a , S, 
with the functional forms: 



(9) 
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-erfc 
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erfc 



-a 2 -5 



aV2 



37re£ft 
32GWe T S 



Ca,8 



2f s 



-erfc 



aV2 

( 3ne 2 ^h \ _ o-2 _ - 
m \ 32Gt^e 2 T: J 2 

ay/2 



( 10 ) 

The transition from the standard regime to this feedback- 
regulation one is shifted toward lower £ for high A4 or, 
in a more modest way, small h. 

3. COMPARISON TO OBSERVED SF LAWS 

3.1. Local spirals and the Small Magellanic Cloud 

In typical local spiral galaxies, large-scale Hi 
and CO r eservoirs have velo city dispersions of 6 — 
10 km s _1 (jCombes et al. 2002) which, for temperatures 
of 10 3 ~ 4 K, corresponds to M « 1 and h « 100 pc. With 
these values, our model naturally reproduces the break 
in SF law observed in entire galaxies and smaller regions 
(but yet > h), as illustrated in Fig. [2] 



erfc : £ i— > 



cxp(-t 2 )dt = l-erf(£). 



Star formation laws and thresholds from ISM structure and turbulence 



3 



\_ h 
! Pa 



-M = 0.6 (a = 0.50) 

-Ji = 1.0 (ct = 0.75) 

-Ji = 1.5 (a = 1.00) 

„« = 3.4 (u = 1.50) 

-Ui = 8.5 (a = 2.00) 



= 100 pc 
= 10 cm" 



/ /// 

//' 
/ / / 

' \ ( \ ( \ i i 



(a) 



■h = 50 pc 
■h = 100 pc 
h = ^00 pc 
h = 500 pc 
■h = 800 pc 



Ji 
Pa 



10 cm 



/ 



/ / / / / 
/ / / / / 
/ / / / / 

A i ^ i i' I 1 1 1 i 1 1 1 1 1 



(b) 



2 
1 



-1 

-2 

-3 
-4 
2 
1 


-1 
-2 

-3 
-4 
2 
1 

-1 
-2 

3 

-4 
2 
1 



-1 

-2 

-3 
-4 

1 2 
log(E [M pc" 2 ]) 

Fig. 1. — Surface density of the star formation rate, computed 
from equation J8), with e = 0.01. The dependences with the Mach 
number M (via the width a of the PDF, panel a), the thickness h 
(panel b) and the density threshold po (panel c) are shown. Panel 
(d) shows the regulation of star formation due to the stellar feed- 
back (solid lines, equation [9} , compared to no regulation, as above 
(dashed lines). The black dashed line is the same in all panels. 
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Fig. 2. — Star formation rate surface density computed with pa- 
rameters (M, h in pc, po m cm -3 ) representative of spiral galax- 
ies (black) and of the Small Magellanic Cloud (SMC, red). The 
models are com pared to observati onal data of the THINGS survey 
(black contours, [Bigicl ct al. 2008), other s piral galaxies (triangles , 
IKennicuttj|19981 ). M51 region s (black dots, Kcnnicut t'et al.ll2007t) 
and the SMC (red contours, IBolatto et al.ll201U ). 



The blue and 



green lines are alternative fits of lower quality (see text). 

The SMC has different properties: we estimate its 
thickness to be ~ 500 pc, but note that its interaction 
with th e Large Magellan ic Cloud could make it even 
thicker (|Besla et al.ll2012j ). Furt hermore, the low m etal- 
licity of the SMC (~ 0.2Z Q , e.g. IBolatto etaI1[2008l) im- 
plies a less efficient cooling of the ISM leading to a tran- 
sition to supersonic turbulence at po = 100 cm" 3 (using 
coolin g calculations similar to those of iBournaud et all 
2010). Using these values, our model matches the obser- 
vations (Fig. [21 red). 

The (relative) uncertainty on h allows us to adjust the 
triplet of parameters to explore other possibilities. First, 
considering that po is independent of the metallicity (and 
thus has the same value as for spirals) would require ad- 
justing the thickness to an unrealistic value (3 kpc) to 
match the data (green curve in Fig. \2§- Second, not ac- 
counting for the structure of the ISM (A4 and h) and only 
focusing on the role of the onset of supersonic turbulence 
through a high po would lead to a too steep relation (blue 
line in Fig. [5]). Alternatively, if we would assume that the 
formation of H2 triggers SF instead of the transition to 
supersonic turbulence, we could model this with an even 
higher threshold (po ~ 100 — 300 cm~ 3 for the forma- 
tion of molecules at subsolar metallicity), and thus get a 
steeper slope failing to fit the observations. Hence, the 
observed SF law in the SMC is best explained by the 
onset of supersonic turbulence and the associated ISM 

4 The thickness of the SMC is derived using lCoinbes et al.l H20021 . 
Chapter 1) with a velocity dispersion of 20 km s —1 and a ro- 
tation speed of 40 km s —1 at the Hi half-mass radius (1.5 kpc, 
IStanimirovic et al]l2004l 1. 
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Fig. 3. — Density PDF of a merger with the parameters adopted 
throughout the paper u\ = <T2, 5 = 3 and m = 0.2 (see text). 

structure^. 

3.2. Application to starbursting mergers 

When they experience a starburst phase, galaxy merg- 
ers often have a part of their gas compressed in a nu- 
clear disk, which we model by converting a mass fraction 
m of the initial PDF into the form of a denser com- 
ponent (using a dimensionless parameter 6, see Fig. [3]): 
(1 — m)f ax (x) + mfa 2 (x/ exp 5). Hydrodynamic simula- 
tions of mergers found such excess of high-density com- 
ponents co mpared to a log- norm al PDF, even out side nu- 
clear disks (jTevssier et al.ll2010HBournaud et al.ll2011aD . 
and PDFs from these works are well fitted by our arbi- 
trary model using S = 3 and m — 0.2. With feedback, 
we obtain 



J SFR 



s 3/2 

[(i 



(1 — m)b ai ^ + mexp 

m)c au0 + mexp (6)c a2t $\ . 



35 

y ]»<t 3 ,s 



(11) 



Such SF law yields a single break, at lower surface den- 
sities than in isolated galaxies, and a flattening at high 
surface densities due to stellar feedback. 

3.3. Clouds, disks and mergers at low and high redshift 
High-redshift (z = 1 — 2) disks have strong turbulence 
(50 - 100 km s _1 , iForster Schreiber et al.1 [2001) lead- 
ing t o a thicker disk (~ 1 kpc. lElmegreen fc Elmcgrccn 
2006), and a higher £ than local spirals. In such a 
regime, feedback induces a transition of the power-law 
index from 3/2 to unity, which best matches the data 
(Fig. 2]). Note that the onset of supersonic turbulence 
probably occurs at a lower density threshold than that 
of nearby spirals (10 cm -3 ). However, this would mainly 
modify the shape of the break, in a range of surface den- 
sities much lower than that of the (non-resolved) high-z 
galaxies considered here. 

The ISM in starbursting mergers is at least twice 
more t urbul e nt than in isolated galax ies (at low- redshift: 
I Irwin! Il994t lElmegreen et al.l 119951 and high-redshift: 

5 In these regimes of low M and Esfr (spiral disks and SMC), 
the feedback has a very mild effect, not affecting the relations. 



Bou rnaud etldl l2011bO . This can increase the scale- 
height of gas reservoirs, although gas in the most central 
regions is not necessarily thicker t han in disk galaxies, 
becau se of its high surface density (jDownes fc Solomon! 
1998). In Fig. 01 we adopt doubled M. and h compared 
to disks at the same redshift. In these dense and turbu- 
lent systems, the regulating role of feedback is already 
significant at z — and increases with z, explaining 
the deviation from the 3/2-slope toward the shallower 
slop es (1.2 - 1.3) as o bserved by iDaddi et all ()2010bf ) 
and lGenzel et al.l (|2010). Furthermore, the same authors 
noted that SF relations globally raised by ~ 0.9 dex for 
mergers compared to spirals, which our model retrieves 
well, thanks to the extended PDF (excess of dense gas 
and strong turbulence). 

The strong turbulence in molecular clouds of the 
Milky-Way (MW) makes them lie close to the low-z 
mergers. Although our model may slightly underpredict 
Ssfr for these objects, its qualitative behavior compared 
to spirals and the general trend due to feedback regula- 
tion are in good agreement with the observations. In fact, 
these molecular clouds are distributed between the model 
predictions with and without feedback, which could indi- 
cate a mix of young clouds not yet affected by feedback, 
and of old ones undergoing strong regulation. 

In summary, the regulation by feedback, the turbu- 
lence increasing with the redshift, the extended PDF 
for starbursting mergers and the higher threshold at low 
metallicities are necessary aspects to retrieve the diver- 
sity of the SF laws. Molecule formation could be con- 
com itant with SF, withou t being the main physical trig- 
ger (|Glover fc Clarkjl20l2l) . 

4. ALTERNATIVE INTERPRETATIONS, POSSIBLE TESTS 

In the proposed framework, a local volume density 
threshold is required to explain an observable surface 
density break. We have proposed that the density thresh- 
old corresponds to the onset of supersonic turbulence, 
generating shocks that trigger the gravitational instabil- 
ities leading to SF. However, other interpretations are 
possible, keeping the same formalism. 

First, on top of self- gravity, the galactic environment 
of a cloud can modify its equilibrium, delaying or acceler- 
ating the onset of SF. The tidal fie ld and the shear tend 
to sta bilize clouds against collapse (|Elmegreen fc Hunterl 
2006). Therefore, the threshold for SF could corre- 
spond to the critical density needed to initiate collapse. 
In pressure equilibrium, p$ would be the density for 
which tidal, centrifugal (shear) and Coriolis forces bal- 
ance self-gravity. Although it depends on the galacto- 
ccntric radius, this tidal density averaged over the galaxy 
is ~ 10 cm -3 , with limited variations with galaxy type 
and redshift. Second, support by magnetic fields might 
stabilize molecular clouds, leading to higher thresholds 
and shallower slopes in all regimes. 

Both these alternative interpretations would translate 
into mild variations of the density threshold with the 
galactic environment. However, if, as we have proposed 
here, the local density threshold corresponds to the on- 
set of supersonic turbulence, highly turbulent media (e.g. 
mergers and high redshift disks) would have very low 
thresholds leading to breaks in the observed SF relation 
shifted to low surface brightness, close to the detection 
limits and thus difficult to detect. Preliminary obser- 
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Fig. 4. — Comparison of our models, using physically motivated parameters, with observations of j 
and high rcdshift (2). 



Jaxies or regions of galaxies at low 



vations of r esolved regions of low -z mergers (Belles et 
al., in prep, iBoauien et all l201l|) started to probe this 
regime (see Fig. [4j : some regions are spiral-like (as ex- 
pected since not all mergers are starbursting) but ~ 20% 
of them have significant Esfr for very low S, well below 
the break of spirals. Larger samples are needed to tell 
apart starbursting regions from spiral-like ones, and to 
confirme the trend suggested here. This could also be 
tested with resolved observations of z — 2 disks, to fur- 
ther probe the physical origin of the local density thresh- 
old. 

5. SUMMARY AND CONCLUSION 

In this Letter, we present an analytical formalism aim- 
ing at describing the relation between the surface density 
of gas and the surface density of SFR observed in several 
types of galaxies. The only two ingredients of the model 
are the gas density PDF (representing the turbulence- 
driven structure of the ISM), and a local SF law, with a 



threshold due to the onset of supersonic turbulence, plus 
regulation by stellar feedback. Our main findings are: 

• When integrated over regions of galaxies or entire 
galaxies, the threshold in the local SF law trans- 
lates into a break at low surface densities. 

• Above this break, the local SF law directly imprints 
the global SFR surface density, leading to a 3/2- 
index power-law, followed by a unity slope at high 
surface density, when feedback takes over turbu- 
lence as the main regulation agent. 

• The slow nature of star formation in nearby spirals 
with Gyr-long gas depletion timescales is mostly 
explained by turbulent regula tion (i.e. diver - 
gent flows disrupting gas clouds, EJmcgrcen 2002). 
Feedback regulation becomes more prominent in 
denser and more turbulent systems (e.g. high- 
z galaxies or Galactic molecular clouds) and pro- 
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duces the observed shallower relation for the entire 
disk population (average slope 1.1-1.3). 

• The different properties of the SMC can explain 
the onset of supersonic turbulence at higher surface 
densities than in spirals, and our model retrieves 
the observed longer depletion times without neces- 
sarily invoking molecule formation or shielding by 
dust. 

• In starbursting mergers, the strong interaction- 
induced turbulence results in high-efficiency SF, 
with a break shifted to low densities but tentatively 
observed in resolved galaxies. 

By simply using the typical values of ISM properties 
from the literature for several galactic environments, our 
model naturally explains the observations, without fit- 
ting them nor fine tuning the parameters. 

The temperature of the gas shielded from ultraviolet 
radiations and in thermal equilibrium would drop when 
the d ensity reaches s» 0.8 cm -3 ([Audit fc Hennebelld 
120101 ) , allowing ISM turbulence to become supersonic at 



these relatively low volume densities. In spiral galax- 
ies, such low threshold would translate into an appar- 
ent break at < 1 M Q pc~ 2 . Instead, according to mod- 
els including not only ultraviolet radiative transfer but 
also turbulent compression, the gas slowly cools below 
~ 10 K and enters the superson ic regime at about 
10 cm -3 (jAudit k Hennebellejl2010f) . which naturally ex- 
plains the observed break at ~ 10 M Q pc~ 2 . 

Our results emphasize that the supersonic turbulence 
plays a key role in triggering and regulating galactic- 
scale star formation, along with feedback in dense and 
turbulent media, and that the resulting ISM structure 
can broadly explain the observed star formation laws and 
thresholds. 
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